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Context. We studied numerically the formation of giant planet (GP) and brown dwarf (BD) embryos in gravitationally 
unstable protostellar disks and compared our findings with directly-imaged, wide-orbit 50 AU) companions known 
to-date. 

Aims. The viability of the disk fragmentation scenario for the formation of wide-orbit companions in protostellar disks 
around (sub-)solar mass stars was investigated. The particular emphasis was paid to the survivability of GP/BD embryos 
formed via disk gravitational fragmentation. 

Methods. We used numerical hydrodynamics simulations of disk formation and evolution with an accurate treatment of 
disk thermodynamics. The use of the thin-disk limit allowed us to probe the long-term evolution of protostellar disks, 
starting from the gravitational collapse of a pre-stellar core and ending in the T Tauri phase after at least 1.0 Myr of 
disk evolution. We focused on models that produced wide-orbit GP/BD embryos, which opened a gap in the disk and 
showed radial migration timescales similar to or longer than the typical disk lifetime. 

Results. While disk fragmentation was seen in the majority of our models, only 6 models out of 60 revealed the 
formation of quasi-stable, wide-orbit GP/BD embryos. The low probability for the fragment survival is caused by 
efficient inward migration/ejection/dispersal mechanisms which operate in the embedded phase of star formation. We 
found that only massive and extended protostellar disks (^0.2 M©), experiencing gravitational fragmentation not only 
in the embedded but also in the T Tauri phases of star formation, can form wide-orbit companions. Disk fragmentation 
produced GP/BD embryos with masses in the 3.5-43 Mj range, covering the whole mass spectrum of directly- imaged, 
wide-orbit companions to (sub-) solar mass stars. On the other hand, our modelling failed to produce embryos on orbital 
distances ^170 AU, whereas several directly-imaged companions were found at smaller orbits down to a few AU. Disk 
fragmentation also failed to produce wide-orbit companions around stars with mass ^0.7 M©, in disagreement with 
observations. 

Conclusions. Disk fragmentation is unlikely to explain the whole observed spectrum of wide-orbit companions to (sub- 
) solar-mass stars and other formation mechanisms, e.g., dynamical scattering of closely-packed companions onto wide 
orbits, should be invoked to account for companions at orbital distance from a few tens to ~ 150 AU and wide-orbit 
companions with masses of the host star < 0.7 Mq. Definite measurements of orbit eccentricities and a wider sample 
of numerical models are needed to distinguish between the formation scenarios of GP/BD on wide orbits. 
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1. Introduction 



With the detection of giant planets (CPs) and brown 
dwarfs (BDs) on orbital distances of the order of tens 



to hundreds AU fe.g. iMarois et all 120081 : iKalasl . I2QQ8I : 

008), 



■ iLafrenierd . 120101 : ISchmidt et al 



"2008), gravitational insta- 
bility and fragmentation of protostellar disks has gained 
a renewed interest as one of the likely mechanisms 
that can explain the fo rmation of wide-orbit co mpanions 
e.g. iBossl. [2001 l201ll: Stamat ellos fc Whitwort h. 2009a; 
Bolev, 2009; Dodson- Robinson et al., 2009; Kratter et alj, 
2010a, b; Vorobyov fc Basu. 2010b). Measurements of disk 
masses in the T Tauri phase of stellar evolution yielded a 
few candidates with disks masses as massive as 0.1-0.5 M^?) 
(lEisner et al.L l2008l : [Andrews etld] . 120091 : llsella et al.L 
1200 9). Such massive disks are expected to be more fre- 
quent in the earlier, embedded phase of star formation (e.g . 
lEisner et al] l2005l : [jQrgensen et al.L 120091 : 1 VorobvQvll2011bl : 
lEisner','2012), reinforcing the viability of disk fragmentation 
at least in this early evolutionary stage. 



Numerical hydrodynamics simulations and semi- 
analytical studies seem to converge on that disk frag- 
mentation is feasible at distances greater than a few 
tens of AU from the c entral star (Maver et al.", "20071; 
Bolevl. l2009URice et all 120101 : |Vorobyov fc Basu, 2010ai 
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Zhu et all . 2012[ ) , where the cooling time becomes s horter 

than the local dy namical timescale ( Gammi^ . l200ll : 

iJohnson fc Gammiel . 120031 : iRice et al.L 120031). Although the 
exact distance beyond which disk fragmentatio n may oper- 
ate is still under deb ate (Rafikoy. 2005 ; Nero fc Bjorkmanl . 
I2009L iMeru fc Bat3 . 1201 iL 1201 2[ ). it now becomes evident 
that parental cores must have enough mass and angular 
momentum to form extended disks with mass sufficient 
to drive the local Toomre Q-parameter be low unity (e.g. 
IVorobvov fc Basul . l2010aL IVorobvovL 1201 Ibl ). 

It has recently become evident that the formation of 
massive fragments via disk gravitational fragmentation 
does not guarantee that the fragments will ultimately 
evolve into GPs/BDs on wide orbits. Gravitational insta- 
bility in the embedded phase of star formation is strong, 
fuelled with a continuing infall of gas from a parent cloud 
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core, and resultant gravitational and tidal torques from spi- 
ral arms are rampant. As a consequence, the majority of the 
fragments are torqued into the inner disk regions, probably 
producing a population of close-in terrest rial and icy giants 
due to tidal disruption/tidal downsizing (iNavakshinl 120101 : 
iBolev et all l201ol : ICha fc NavakshinL l201lD, or driven di- 
rectly onto the star, triggering intense accretion and lumi- 
nosity bursts similar in magnitu de to FU Orionis eruptions 
(|Vorobvov fc Basul , l2006l l2010al ) . Some of the fragments are 
dispersed by tidal torques exerted by the spiral arms before 
they can dissociate molecular hydroge n in their in t eriors 
and cont ract to planetary-si zed ob j ects (iBolev et aD , I2OIOI : 
IVorobvovi , i2011ai : iNav a kshin et al.L l201l[ ). A few fragments 
may be scattered away from the disk via many-body grav- 
itational interaction with outer fragments or fully formed 
sub-st ellar objects, producing fre ely floa ting sub-stellar ob- 
jects (jStamatellos fc Whitworthl . [2009al: Ib asu fc Vorobvovl 
I2OI2D . 

Therefore, the question of whether fragments can settle 
into stable orbits at distances where the fo rm, fro m several 
tens t o hundreds AU, remains to be open. iBaru teau et al. 
(|201lf ) studied planet migration in graviturbulent disks 
with mass 0.4 Mq and argued that Jupiter-mass plan- 
ets (or higher) initially placed at 100 AU migrate inward 
on timescales of 10^ yr and are unlikely to stay on wide 
orbits. Michael et aJ] (|201ll ) found even faster migration 
timescales of 10^ yr for a Jupiter-mass planet in a 0.14 Mq 
disk, though the planet may stall near the inner Lindblad 
resonance of the dominant spir al mode. Other numeri- 



cal hydrodynamics simu lations (Vorob vov fc Ba su. 2006, 
l2010ai : ICha fc Navakshinl . I2OIII : iMachida et al.L 12011) also 



revealed fast inward migration of the forming fr a gments 
in the disk. On the other hand, I Vorobvov fc Basul (l2010bD 
studied the long-term {r^ several Myr) evolution of frag- 
menting protostellar disks in the thin-disk limit and found 
that while most disks indeed fail to produce stable compan- 
ions on wide orbits, in agreement with previous studies, a 
small subset of models can form GPs at distances of the 
order of tens to hundreds AU. They concluded that only 
those fragments that happen to form in the late embedded 
phase, when gravitational instability and associated torques 
are getting weaker, may open a gap in the disk and mature 
into GPs on wide orbits. The low probability for survival 
of the fragments formed vi a disk grav i tatio nal fragmenta- 
tion was also confirmed by IZhu et aH (j2012l ) , who studied 
numerically two-dimensional disks subject to mass loading 
and found that only 3 fragments out of 13 became massive 
enough to open a gap in the disk and essentially stopped 
migrating. 

In this pa p er, w e improve the model of 
I Vorobvov fc Basul (|2010b ^ by including a detailed ther- 
mal balance in the protostellar disk, t hus remov ing the 
barotropic relation closure adopted in I Vorobvov fc Basul 
(|2010b[ ) which is kno w n to o yerpredict the number of 
fragme nts (e.g. iBatd . 120091 : IStamatellos fc Whitworthl . 
l2009bD . In our numerical hydrodynamics simulations, 
we form disks self-consistently during the gravitational 
collapse of pre-stellar cores and not introduce them 
artificially. This allows us to determine the range of masses 
and angular momenta in prestellar cores for which disk 
fragmentation and formation of wide-orbit companions 
can take place. We also avoid replacing fragments with 
point sink particles, thus studying the evolution of GP/BD 
embryos rather than fully formed planetary- or sub-stellar 



objects. This is done to avoid a premature introduction 
of sink particles, the effect that can influence the number 
of surviving fragments due to essential indestructibility of 
point-sized objects. On the other hand, the formation of 
planetary /sub-stellar-sized objects is not resolved in the 
current approach, which may affect the short-term surviv- 
ability of the fragments. We focus on models that produce 
planetary- or sub-stellar-mass embryos on quasi-stable 
orbits at radial distances where disk fragmentation takes 
place (> 50 AU). 

The organization of this paper is as follows. A brief de- 
scription of the numerical model is provided in Section |2j 
The parameter-space study of disk fragmentation and frag- 
ment survival is provided in Section |3l Formation of GP 
and BD embryos in wide orbits is presented in Sections |4] 
and ini respectively. Feasibility for the formation of mul- 
tiple companions in wide orbits is discussed in Section |6l 
Comparison of numerical results with observed sub-stellar 
objects in wide orbits is performed in Section [71 The main 
results are summarized in Section [51 

2. Model description 

Our numerical model is described in detail in 
IVorobvov fc Basul (|2010b ^ and is briefly reviewed be- 
low for the reader's convenience. We use numerical 
hydrodynamics simulations in the thin-disk approximation 
to compute the gravitational collapse of pre-stellar cores of 
various initial mass and angular momentum. This approx- 
imation is an excellent means to calculate the evolution 
for many orbital periods and m any model parame t ers, an d 
its justification is provided in ([Vorobvov fc Basul . l2010al ). 
To avoid too small time steps, we introduce a "sink cell" 
at Tsc = 6 AU and impose a free boundary condition such 
that the matter is allowed to flow out of the computational 
domain but is prevented from flowing in. The sink cell 
is dynamically inactive; it contributes only to the total 
gravitational potential and secures a smooth behaviour of 
the gravity force down to the stellar surface. We monitor 
the gas surface density in the sink cell and when its value 
exceeds a critical one for the transition from isothermal 
to adiabatic evolution, we introduce a central point-mass 
object. 

The simulations continue into the embedded phase of 
star formation, during which a protostellar disk is formed. 
In this stage, the disk is subject to intense mass loading 
from the remnant of the initial pre-stellar core — the so- 
called envelope. The self-consistent disk-envelope interac- 
tion is a key feature of our model allowing us to observe 
repetitive episodes of disk fragmentation in some models. 
During the diks evolution, 90% of the gas that crosses the 
inner boundary is assumed to land onto the central ob- 
ject plus the sink cell. The other 10% of the accreted gas 
is assumed to be carried away with protostellar jets. The 
simulations are terminated in the late T Tauri phase after 
more than one Myr of disk evolution when nearly all of the 
envelope material has accreted onto the resulting star plus 
disk system. 

Models presented in this paper are run on a polar coor- 
dinate (r, (/)) grid with 512 x 512 zones. The radial points 
are logarithmically spaced, with the innermost cell outside 
the central sink having size 0.07-0.1 AU depending on the 
cloud core size (i.e., the radius of the computational re- 
gion). The latter varies in the 0.025-0.12 pc (5000-24,000 
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AU) limits. The radial and azimuthal resolution are ^1.0 
AU at a radial distance r ^ 100 AU. 

The Truelove criterion states that the local Jeans length 
must be resolved by at least four numer ical cells in or- 
der to correctly capture disk fragmentation ([Truelove et al.L 
I1999D . For a circular fragment of radius Rq with surface 
density inversely proportional to radius, H = T^oRo/r^ the 
kinetic energy due to random motions can be expressed as 



(1) 



where is the surface density at the fragment-disk inter- 
face. The velocity dispersion of a thin disk with two trans- 
lational degrees of freedond is {v'^) = 27^To//i, where IZ is 
the universal gas constant, /i = 2.33 is the mean molecu- 
lar weight, and Tq is the gas midplane temperature in the 
fragment. 

The corresponding gravitational energy of the fragment 

is 



/■ 



= —27T rg^T^rdr = 



-27r^GI]Qi?Q, 



(2) 



where we have taken in to account that Qr = ttGE for a 
disk with S oc r"^ fsee iBinnev fc Tremain^ . I1987L p. 77). 
The resulting Jeans length i?j is calculated from the virial 
theorem 2£^kin + = as 



ttGSo' 



(3) 



Fragments usually condense out of densest sections of spiral 
arms. The typical surface densities and temperatures in spi- 
ral arms do not exceed 100 g cm~^ and 10 K (see Figures [2l 
[H [8]in this paper and IVorobvovl (|2011bf )). Adopting these 
values for Sq and Tq, the corresponding Jeans length is 
i?j ^ 20 AU. 

In models showing disk fragmentation, the radial and 
azimuthal grid resolution at r = 100 AU is ~ 1.0 AU and 
the Jeans length is resolved by roughly 20 grid zones in 
each coordinate direction. On our logarithmically spaced 
radial grid, the Truelove criterion is expected to break only 
at r ^ 500 AU where the grid resolution starts to exceed 
5.0 AU. Fragmentation takes place mostly at radial dis- 
tances from a few tens to a few hundreds AU. Fragments 
that are seen in our models at larger distances are most 
likely scattered from the inner disk regions due to gravita- 
tional interaction with other fragment^. The radii of the 
survived fragments (see the ninth column in Table [2]) lie 
between 10 and 20 AU, implying that the fragments are re- 
solved on the two-dimensional mesh by at least 30-60 grid 
zones in the inner 500 AU. We therefore conclude that the 
numerical resolution in our models is sufficient to capture 
disk fragmentation correctly. On the other hand, contrac- 
tion of the survived fragments to planetary-sized objects 
cannot be modelled in the current approach. This may 
have consequences for the survivability of the fragments. 



^ Our disks are not razor thin, but are characterized by a 
vertical scale height increasing outward, as occurs in flared disks. 
Taking into account the vertical degree of freedom changes the 
corresponding Jeans length by a factor of only 1.5. 

^ In the most extreme cases, some of the fragments may 
even be ejected from th e disk into the intracluster medium 
(|Basu Vorobvovi [201^ 1. 



resulting in an increased probability of tidal destruction 
of AU-sized objects as compared to planetary-sized ones. 
Accreting sink particles are needed to correctly follow the 
evolution of fully formed giant planets and brown dwarfs. 

2.1. Basic equations 

In the paper of IVorobvov fc Basul (|2010a[ ), a barotropic 
equation of state was used to close the equations of hydro- 
dynamics. In this work, we include detailed thermal physics 
in our model, the main concepts of which are briefly re- 
viewed below. The basic equations of mass, momentum, 
and energy transport are 



(v-n)^. 



(4) 



(5) 



as 
dt 

d_ 
dt 

de 

g^^'^p- {evp) = -ViSp • i;^) - A + F + (Vi;)^^, : H^^. , (6) 

where subscripts p and p' refers to the planar components 
(r, (j)) in polar coordinates, S is the mass surface density, e 
is the internal energy per surface area, V is the vertically 
integrated gas pressure calculated via the ideal equation 
of state as 7^ = (7 — l)e with 7 = 7/5, Z is the radially 
and azimuthally varying vertical scale height determined 
in each computational cell using an assumption of local hy- 
drostatic equilibrium, Vp = VrV -\- V(f)<j) is the velocity in the 
disk plane, Qp = grV -\- g(/y(i) is the gravitational accelera- 
tion in the disk plane, and Vp = rd/dr -\- 4>r~^d/d(j) is the 
gradient along the planar coordinates of the disk. We note 
that the adopted value of 7 neglects a possible stiffening of 
the equation of state at low temperatures (< 100 K) where 
the ratio of specific heats may approach a value typical fo r 
a monatomic gas, 7 = 5/3 (Masunaga fc Inutsukal , l2000f ). 
The adopted equation o f state may be ini portant for disk 
gravitational instability (iBoley et al.L |2007[ ) , though the re- 
cent study of lZhu et al.l (|2012[ ) found that fragmentation of 
two-dimensional disks does not depend sensitively on 7 in 
the range from 7/5 to 5/3. 

Turbulent viscosity is taken into account via the viscous 
stress tensor 11 expressed as 



U = 2T.v [Vv - -(V • i;)e 



(7) 



where e is the unit tensor and Vv is the symmetrized 
velocity gradient tensor. We parameterize the magnitude 
of kinematic viscosity v using a modified form of the a- 
prescription 



(8) 



where Cg = jT, is the square of effective sound speed cal- 
culated at each time step from the model's known V and 
E. The function Ta,{T) = 27r~^ tan""*^ [(^d/^)"^^] is a modifi- 
cation to the usual a-prescription that guarantees that the 
turbulent viscosity operates only in the disk and quickly 
reduces to zero beyond the disk radius r^- In this paper, we 
use a spatially and temporally uniform a, with its value set 
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to 5 X 10~^ to take into account mass and angular momen- 
tum transport via such mechanisms as the magnetorota- 
tional instability. Transport of mass and angular momen- 
tum via gravitational instability is self-consistently taken 
into account via solution of the Poisson equation for the 
gravitational potential of the disk and envelope. 

The radiative cooling A in equation (|6]) is determined 
using the diffusion approximation of the vertical radiation 
transport in a one-zone mod el of the vertical disk structure 
(jJohnson fc Gammid . [20031 ) 



A = J^eCrT^ 



(9) 



where <j is the Stefan-Boltzmann constant, T is the mid- 
plane temperature of gas, and = 2 + 20 tan~^ (r) / (Stt) is 
a function that secures a correct transition between the op- 
tically thick and optically t hin regimes. We use frequency- 
integrated opacities of iBell fc LinI (|199T| ). The heating func- 
tion is expressed as 



r = J^nCr T-tr 



(10) 



where Ti^r is the irradiation temperature at the disk sur- 
face determined by the stellar and background black-body 
irradiation as 



i^irr(r) 



(11) 



where Tbg is the uniform background temperature (in our 
model set to the initial temperature of the natal cloud core 
^init = 10 K) and Firr(r) is the radiation flux (energy per 
unit time per unit surface area) absorbed by the disk sur- 
face at radial distance r from the central star. The latter 
quantity is calculated as 



i^irr(r) 



cos7ij 



(12) 



where 7irr is the incidence angle of radiation arriving at the 
disk surface at radial distance r. 

The stellar luminosity is the sum of the accretion 
luminosity L^^acr = GM^M /2R^^ arising from the gravi- 
tational energy of accreted gas, and the photospheric lu- 
minosity L^^ph due to gravitational compression and deu- 
terium burning in the stellar interior. The stellar mass 
and accretion rate onto the star M are determined 
self-consistently during numerical simulations using the 
amount of gas passing through the sink cell. The stel- 
lar radi us is calcula t ed us in^ an approximation for- 
mula of iPalla fc Stahlerl (|l99lf ). modified to take into ac- 
count the formation of the first molecular core. The pho- 
tospheric luminosity I/*^ph is taken from the pre-main se- 
quence tracks for the l o w- ma ss stars and BDs calculated by 



iD^Antona fc Mazitelhl (119971). More details on the numeri- 
cal code are given in IVorobvov 



Basul (|2010aD . 



2.2. Initial conditions in pre-stellar cores 

We considered two limiting cases to describe the initial dis- 
tribution of the gas surface density H and angular velocity 
Tl in the pre -stellar cores. The first distribution, taken from 
iBasul (|1997[ ). is typical of pre-stellar cores formed as a re- 
sult of the slow expulsion of magnetic field due to ambipolar 



diffusion, with the angular momentum remaining constant 
during axially-symmetric core compression 



S = 



__ro^o_ 




(13) 



(14) 



Here, Qq and are the angular velocity and gas surface 
density at the disk center and ro = ^/Ac^/irGT^o is the ra- 
dius of the central plateau, where Cg is the initial sound 
speed in the core. The gas surface density distribution de- 
scribed by equation ([13]) can be obtained (to within a factor 
of unity) by integrating the three-dimensional gas density 
distribution characteristic of Bonnor-Ebert s pheres with a 
positive density-perturbation amplitude A (jPapp fc Basul . 
2009). In all models the value of A is set to 1.2, except for 
model 4 for which A=3.3. 

The second set of initial conditions described by 
r^=const and I]= const represent the other limiting case sug- 
gested in iBoss fc HartmannI (|200lf ) . Cores that can be de- 
scribed, to a first degree of accuracy, by spatially constant S 
and Q can form via gravitational fragmentation of filamen- 
tary structures, which are often encountered in numerical 
hydrodynamics simulations of the turbulent fragmentation 
of giant molecular clouds. Another possible mechanism for 
the formation of such cores may be the planar compression 
of pre-stellar condensations by shocks and UV radiation of 
massive stars. 



2.3. Tracking the fragments 

Motivated by the absence of sink particles in our grid- 
based code, we have designed a fragment-tracking algo- 
rithm which allows us to follow the trajectory of the frag- 
ments and calculate their physical parameters. The most 
straightforward way to identify fragments in the disk is to 
set a threshold gas surface density Scrit which would help to 
distinguish between the fragments and the rest of the disk. 
We however quickly found out that setting a single value of 
Scrit, independent of radial distance, did not work because 
too small a value for Scrit might result in spiral arms in the 
inner disk being identified as fragments. We therefore need 
to define a radially varying value of Scrit which is greater 
at smaller radii and vice versa. 

Noting that the fragments are usually characterized 
by peak surface densities that are greate r than the local 
azimut hally averaged surface density S (|Vorobvovi l2010l 
l2011b[ ). the threshold gas surface density can be defined 
as Scrit = C'S, where C is an empirically determined con- 
stant. Numerical simulations show that the azimuthally av- 
eraged gas surface density_in gravitationally uns table disks 
declines with radius as S = i;o(ro/r)^-^ (e.g. IVorobvovi 
l2010l : [Rice et al.Ll2010f ). where is the typical azimuthally 
averaged density at tq. For gravitation ally unstable disks , 
So = 20 - 50 g cm-2 at r ^ 100 AU (|Vorobvovl l2011b[ ). 
We therefore have chosen the following expression for the 
threshold density 



^100 



100 AU ^ 



1.5 



(15) 
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The best value for threshold surface density at r = 100 AU 
is found by trial and error method to be Hioo = 200 g cm~^. 
Note that this value is greater than the mean gas surface 
density So by at least a factor of several. The local maxima 
in the gas surface density that exceed Scrit usually represent 
the true fragments rather than the local maxima in the 
disk and/or spiral arms. Although we may occasionally miss 
some of low-density fragments, this does not invalidate our 
main conclusions. 

After the radial Tc and angular coordinates of the lo- 
cal maximum representing the center of the fragment have 
been identified on the computational mesh, we determine 
the neighbouring cells that belong to the fragment by im- 
posing the following two conditions on the gas pressure V 
and gravitational potential ^ 



dV 
dr' 

dr' r' 



r' d^' ^ ' 
1 5£ 



(16) 
(17) 



where r' = r — Vc and = (j) — (j)^. The first condition 
mandates that the fragment must be pressure supported, 
with a negative pressure gradient with respect to the center 
of the fragment. The second condition requires that the 
fragment is kept together by gravity, with the potential well 
being deepest at the center of the fragment. A substantial 
support against gravity may be provided by rotation but 
we assume that this does not invalidate our criteria, i.e., no 
fragments assume a torus shape. 

In practice, we start from the grid cell corresponding to 
the center of the fragment and proceed in eight directions 
(along the four coordinate directions and also at median 
angles to them) until at least one of the above mentioned 
criteria is violated in every direction. This procedure helps 
to identify an approximate shape of the fragment. We then 
check all the remaining grid cells that are encompassed by 
this octahedral shape and retain only those that meet both 
criteria (flB]) and ([17]). In addition, we filter out cells with 
the gas surface density lower than that defined by equa- 
tion (p!5|) . even if these cells still fulfil both criteria. We 
found that such cells are hkely to belong to the circum- 
fragment disk rather than to the fragment itself. The cells 
that belong to the fragment are later utilized to calcu- 
late the mass and Hill radius and also the gravitational 
torque exerted on the fragment, while the cell correspond- 
ing to the center of the fragment is used to calculate the 
trajectory of the fragment. The characteristics of the frag- 
ments thus found depend somewhat on the adopted value of 
Sioo = 200 g cm~^. However, this dependence is not criti- 
cal and, for instance, the estimated masses of the fragments 
change by only about 10% if Sioo is varied by a factor of 
two. 



3. Fragmentation and survival of fragments 

We have considered the evolution of 60 initial pre-stellar 
cores with masses ranging from 0.1 Mq to 2.0 Mq and 
ratios of rotational to gravitational energy j3 lying between 
0.27% and 2.2%. We have also varies the magnitude of the 
initial positive density perturbation A in the 1.2-3.3 limits 
and considered cores with distinct initial radial profiles of 
E and ^. 
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Fig. 1. Phase space of /3, the initial ratio of rotational to 
gravitational energy, vs. core mass Mc. The solid/open cir- 
cles correspond to models with/without disk fragmenta- 
tion. The black region is the region where the formation of 
sub-stellar- or planetary-mass companions on wide orbits 
is found in this study. The dark-shaded region is the re- 
gion in which both fragmentation and ejection events may 
occur, and the light-shaded region is the one where only 
fragmentation events occur. The white area marks the re- 
gion where no fragmentation is observed. Arrows illustrate 
uncertainties associated with a coarse grid of models and 
indicate models in the p : phase space that might have 
shown (with a 50% probability) disk fragmentation (ejec- 
tion/companion formation) had we considered a finer grid 
of models. Each pair of data in parenthesis indicate the 
mean disk radius and mass for the corresponding model. 



For each model we ran our fragment-tracking algorithm 
at various evolution times in order to identify models that 
experienced disk fragmentation. We found that most of 
the models showed disk fragmentation but at the same 
time failed to produce wide-orbit companions. The major- 
ity of the fragments were torqued into the inner disk region 
and through the sink cell (6 AU) producing mass accre- 
tion and luminosity bursts similar in magnitude to those 
of FU-Ori-type objects (Vorobyov & Basi], I2006L l2010aD . 
while the remaining few were ejected from the disk into 
the intrac luster medium via many -body gravitational in- 
teraction j|Bas^^^mbvo^2^2)_or dispersed via tidal 
torques (Vor obvovl l2011al: iBolev et al.L 120101) . IBolev et al.l 
(2010) and .Zhu et al.l (|2012f ) also found support for clump- 
driven FUor events. Only 6 models out of 60 have shown the 
formation of stable companions in wide orbits, with their 
mass ranging from 3.5 to 43 Jovian masses. Ou r ob tained 
survival probability is even lower than that of IZhu et all 
(|2012,) (3 out of 13), though they followed the evolution of 
the fragments for a significantly shorter time period. 

Figure [1] illustrates our findings on the phase space of {3 
versus Mc covered in our modeling. The light-shaded area 
defines the region where disk fragmentation takes place, 
while the dark-shaded area outlines the region where both 
fragmentation and ejection of the fragments may occur. 
The data related to disk fragmentation and ejection of the 
fragments were taken from numeric al hydrodynamics sim- 
ulations of iBasu fc Vorobvovl (|2012t ). In the present study 
we added several models that revealed the survival of the 
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fragraents. The total nuraber of raodels used in the current 
study amounts to 60 but only 16 key models (those lying at 
the boundary of the considered regions) have been shown 
in Figure [1] to avoid overcrowding. 

Previous numerical and analytical studies of protostellar 
disk evolution demonstrated that gravitational fragmenta- 
tion is possible at distances ^ 50 AU ( e.^. Cla rkd 12009 : 



Stamatellos fc Whitworthl . 12009a] : iBolevl . .2009: Rice et al 
201 oh and in disks with mass ^0.1 Mq (e.g. iRice et al 



2003'; 'Mav er et all 120071 ). Disk masses and radii in our 
models can be calculated by disentangling the disk and in- 
falling core on the computational mesh. This is however 
not an easy task. We do this in practice by first construct- 
ing the azimuthally averaged gas surface_density profiles S 
and then applying a threshold value of Sd2e = 0.5 g cm~^ 
between the disk and the core. We also use the radial gas ve- 
locity profile to see where the infalling envelope lands onto 
the disk fsee IVorobvovi [2Qllbl for details). The disk radius 
estimates are further complicated by substantial outward 
excursions (or even complete ejections) of some fragments. 
It is not always clear if the fragment belongs to the circum- 
stellar disk or it has already detached from it and should 
be treated as an external companion (see e.g. Figure [2j). We 
assume that the fragment belongs to the disk if S between 
the fragment and the disk does not drop below 0.05 g cm~^, 
a value typical for the interface between the disk and the 
external environment (|Vorobvovi l2011b[ ). In the opposite 
case, the fragment is supposed to have detached from the 
natal disk and is treated as a separate entity. 

Our choice of a threshold of 0.5 g cm~^ is motivated by 
the fact that the resulting distribution of disk radii peaks 
between 10^ and 10^ AU; adopting an order of magnitude 
lower threshold would shift the entire distribution up by 
about a factor of 2-3 towards bigger values. Sizes of em- 
bedded disks are very poorly constrained by observations 
and are typically assumed to be on the order of 100 AU 
or less based on simple centrifugal radius arguments. In re- 
ality, however, gravitational and viscous transport of mass 
and angular momentum will cause disks to spread to sizes 
greater than the corresponding centrifugal radii. We note 
that there are some Hmited observations sup porting the ex- 
istence of lar^e protostellar disks (e.g. Enoch et al.L 120091 
iJorgensen et al.LT2009f ) . Ultimately, the correct threshold to 
adopt in order to disentangle the disk and core in the sim- 
ulations will remain uncertain until the masses and sizes of 
protostellar disks are better constrained from observations. 

Each pair of data in Figure [1] indicates the mean disk ra- 
dius and mass for the corresponding model marked with cir- 
cles. The mean values are calculated by time- averaging the 
instantaneous values over the duration of the Class I phase, 
in which disk fragmentation is usually most vigorous. It is 
evident that the disk mass and radius both have to exceed 
some threshold values in order for the disk to fragment. In 
particular, models with higher /3 experience disk fragmen- 
tation at lower disk masses then their low-/3 counterparts. 
The minimum mean disk mass at which fragmentation can 

fr 

take place according to our modeling is = 0.07 Mq for 
P ^ 1.3%. The latter value may increase by about a factor 
of 2 for models with lower /3. We take into account here the 
uncertainty illustrated by the errors in Figure [T] associated 

fr 

with a coarse grid of models. Our found values of are 
in re asonable agreement with previous es timates, 0.1 Mq 
(e.g. iRice et all . 120031 : iMaver et all I2007D . The minimum 



radius of a disk experiencing fragmentation in our models 

fr 

is = 140 AU for /3 ^ 0.3%. The latter value may increase 
by up to a factor of 2 for models with higher p. Our derived 

— fr 

values of are difficult to compare with other studies of 
disk fragmentation because the latter tend to provide the 
minimum radial distance at which fragmentation can take 

place (rfr) and not the minimum disk radius (i^^) ^^^^ 

— fr 

required for disk fragmentation to occur. Obviously, R^ 
must be grater than rfr and our simulations suggest that 

r'^ = (3 - 6)rfr. 

The black area in Figure [T] marks the region where the 
formation of GP and/or BD companions on wide orbits is 
found in the present work. None of our models have shown 
the formation of wide-orbit companions for /3 ^ 1.5% but 
we attribute this to a rather narrow sampling of models 
at this region. Therefore, we extrapolated the companion- 
forming domain to /3 ^ 1.5% assuming the same likelihood 
for the formation of wide-orbit companions as for models 
with lower p. The companion- formation domain is notably 
narrower than both the fragmentation and fragmentation 
plus ejection domains, particularly as far as the initial core 
mass is concerned. It appears that Mc has to be greater 
than 1.2 Mq and the ratio of rotational to gravitational 
energy has to exceed 0.5% in order to make the formation 
of wide-orbit companions possible. Moreover, the minimum 
mean disk mass at which the formation of wide-orbit com- 
panions is found in our modeling is M^ = 0.21 Mq for 

models with P > 0.8%. The value of M^' " may increase 
somewhat for lower- /3 models. The minimum mean disk ra- 
dius that is required for the formation of wide-orbit com- 

c.f . 

panions is R^ ' = 370 AU and this value increases with 
increasing p. Our found values of the minimum disk ra- 
dius and mass both depend on the adopted threshold. For 
instance, adopting a factor of five lower density threshold 
Sd2e = 0.1 g cm~^ would yield a factor of 1.2 greater mini- 
mum disk mass M^-^- = 0.25 Mq and a factor of 1.3 greater 

— c.f. 

minimum disk radius R^ = 470 AU. 

It is not clear why models with Mc < 1.2 Mq fail to 
produce wide-orbit companions. Large masses and radii of 
protostellar disks in companion-forming models must cer- 
tainly be factors that help wide-orbit companions to sur- 
vive. Another likely reason why low-Mc models fail to form 
wide-orbit companions is that disk fragmentation in these 
models is mostly confined to the embedded phase of star 
formation and is sustained by continuing mass loading from 
an infalling envelope. A fragment may escape inward mi- 
gration if the net torque exerted on it is not negative. 

To a first order of accuracy, the net torque acting on the 
fragment in a gravitationally unstable disk can be expressed 
as the sum of the gravitational torques from the inner Tin 
and outer T^ut parts of the disk with respect to the current 
position of the fragment. Spiral arms and other fragments 
are the main contributors to the total torque. Due to the 
trailing nature of the spiral arms in gravitationally unstable 
protostellar disks. Tin is usually positive and T^ut is usually 
negative. The condition for the fragment to avoid inward 
migration can then be written as 



dL 
~dt 



Tin + T^ut ^ 0, 



(18) 
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where L is the angular momentum of the fragment. 
Evidently, the fragment will stay in the disk for as long 
as iT^utl < Tin. 

Even if the fragment forms near the disk outer edge 
where inequality (p!8|) is likely to be fulfilled initially, it 
can start migrating inward during the subsequent evolu- 
tion due to continuing mass loading from the parental core. 
The mass infall onto the disk is a double-edged sword effect: 
it pr omotes disk fragmentation by increasing the disk mass 
(e.g. IVorobvov fc Basul . I2QQ6L l2QlQal : iKratter et all . l2QlQaD 
but it also deposits a sub-Keplerian mate r ial at /near the 
disk outer region (e.g. lVisser fc DullemondlJ2QlQf ). The lat- 
ter acts to increase |7^ut| as the mass near the disk outer 
edge accumulates so that the fragment starts to migrate 
inward when |7^ut| becomes greater than Tin- In addition, 
a sub-Keplerian material falling onto the disk can exert a 
torque onto the fragment and push the later in towards the 
star. 

On the contrary, models with Mc ^ 1.2 Mq are usu- 
ally characterized by disks that are sufficiently massive and 
large to experience fragmentation not only in the embedded 
phase and but also in the T Tauri phase when mass load- 
ing onto the disk diminishes (see Figures [2] and [8]) . When 
formed at the disk outer regions where inequality (p!8]) is 
fulfilled, fragments in these models have chances to open a 
gap in the disk and settle on quasi-stable, wide orbits. The 
main conclusion drawn from our parameter space study is 
that disk fragmentation is not sufficient to guarantee the 
formation of GP or BD companions on wide orbits. 

In this paper, we used a spatially and temporally uni- 
form a = 5 x 1Q~^. This choice is based on the work of 
IVorobvov fc Basul (|2QQ9al lbl). who found that models with 
a = 10~^ reproduce well the slope of the mass accretion 
rate-stellar mass relation for young brown dwarfs and low- 
mass stars (~ 1.0 Myr old), though slightly overpredicting 
the mean accretion rates. Greater value of a act to weaken 
the strength of gravitational instability due to an overall 
increase in the efficiency of mass transport and the corre- 
sponding decrease in the total disk mass. However, large 
values of a (^ 10~^) destroy circumstellar disks during less 
than 1.0 Myr of evolution and are thus inconsistent with 
mean disk lifetimes on the order of 2-3 Myr. We therefore 
conclude that a moderate increase in the adopted value of a 
(by a factor of several) can shift the fragmentation bound- 
ary in Figure [1] towards higher values of Mc and /3 but is 
not expected to shut off disk fragmentation completely. 

Below, we present six models that have shown the for- 
mation of GP/BD embryos on wide, quasi-stable orbits via 
disk gravitational fragmentation. We pay most attention 
to models 1-3 and provide only main results for the other 
three models. The parameters of the models are listed in 
Table m 

4. Formation of a planetary-mass companion 

In this section, we describe the formation of an 11- Jupiter- 
mass companion around a 1.2 Mq star. Figure [2] presents 
a series of the gas surface density images showing the evo- 
lution of the disk in model 1 starting soon after the forma- 
tion of the central object {t = 0.05 Myr) and ending in the 
T Tauri phase {t = 2.37 Myr). The box size is 3000 AU 
on each side and represents a small subregion of the over- 
all computational domain. The time elapsed after the for- 
mation the central protostar is shown in each image and 



Table 1. Model parameters 


model 


Mc 




So 












(Mo) 


(%) 


(g cm~^ 


) 






(AU) 


1 


1.7 


0.56 


3.3 X 10" 


-2 


2.4 X 10" 


-14 


3770 


2t 


1.2 


0.88 


1.3 X 10" 


-2 


1.3 X 10" 


-14 




3 


1.5 


0.56 


3.1 X 10" 


-2 


2.6 X 10" 


-14 


3430 


4 


1.4 


0.56 


1.1 X 10" 


-1 


6.0 X 10" 


-14 


1890 


5 


1.55 


1.27 


3.6 X 10" 


-2 


4.0 X 10" 


-14 


3430 


6 


1.4 


0.56 


4.0 X 10" 


-2 


3.9 X 10" 


-14 


3090 



' Models with spatially constant gas surface density S = So and 
angular velocity Q = Qo- 

3 2.4 1.8 1.2 0.6 O -0.6 -1.2 







1000-^ 


0.05 Myr 




0.1 Myr 


0.22 Myr 


-1000^ 




•4 






1000-^ 

s \ 

^-1000^ 

1 


0.3 Myr 


0.43 Myr 


0.64 Myr 

• 


Radial distant 


0.93 Myr 


1.4 Myr ^ 


1.8 Myr 


1000-^ 


2.1 Myr 




2.18 Myr 


2.37 Myr 












-1000^ 











-1000 1000 -1000 1000 -1000 1000 



Radial distance (AU) 

Fig. 2. Gas surface density distribution in model 1 shown at 
various times since the formation of the central protostar. 
Only the inner 3000 x 3000 AU box is shown, the total 
computational region extends to 22000 AU. The scale bar 
is in log g cm~^. Note a fragment on a stable orbit in the 
bottom row. 



the minimum gas surface density plotted in the Figure is 
0.06 g cm~^ or -1.2 in the log units. 

The forming disk is gravitationally unstable and first 
fragments start to appear in the disk as early as at 50 kyr 
after the formation of the central star. Gravitational per- 
turbations from spiral arms and massive fragments cause 
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1.0 

Time (Myr) 

Fig. 3. Number of fragments vs. time in model 1. The num- 
ber of fragments at a given time instant is calculated using 
the fragment tracking algorithm described in § 12.31 An in- 
crease in the number of fragments shows recent fragmenta- 
tion, and a decrease shows recent destruction/accretion of 
the fragments. 



significant radial motions in the disk. As a result, the 
disk appears as a somewhat chaotic structure with dense 
filaments connecting the fragments. Nevertheless, a near 
Keplerian rotation can be retrieved after azimuthal aver- 
aging. Fragmentation is predominantly concentrated to the 
intermediate and outer disk regions, a consequence of the 
mass infall and stellar irradiation, and no fragmentation is 
evident at r < 50 - 100 AU. 

One may notice from Figure [2] that the number of 
fragments is varying with time, indicating that the frag- 
ments may be tidally destroyed or otherwise lost by the 
disk. For instance, fra gments may mig rate inward onto 
the star (Vorobyov & Basu, 2010a; Machida et al.', '2011) 
or get ejected from the d isk into the intracluster medium 
(|Basu fc Vorobvovl l2012f ). As a result of these migra- 
tion/ejection/destruction processes, only one fragment sur- 
vives after 1.7 Myr of the disk evolution. The bottom row 
in Figure [2] reveals the typical picture with the surviving 
fragment opening a gap and inducing spiral waves in the 
disk. The fragment itself is connected with the inner and 
outer disks by a wake of enhanced surface density which 
trails/leads the fragment outside/inside its orbit. 

Figure [3] presents the number of fragments in the disk 
at a given time instant. The number of fragments varies 
with time and a maximum value {N^^^ = 7) is reached at 
t ^ 0.2 Myr and t ~ 0.9 Myr. Between these two maxima, 
a local minimum with just two fragments in the disk occurs 
at t ~ 0.5 Myr. The disk mass at the end of the embed- 
ded phase (t=0.65 Myr) is ^ 0.3 M©, sufficient to sustain 
fragmentation in the disk for another 0.8 Myr. 

Figure [4] shows four zoomed- in images of the surviving 
fragment taken during a time period of 1180 yr covering 
just two orbital periods of the fragment. The box size is 
500 X 500 AU. The fragment is outlined by the yellow curve, 
which is found using the tracking algorithm described in 



3 
< 

o 
c 

B 

I 400 
CO 

300 





2.27 Myr 


^^^^^^^^ 


2.2723 Myr 






2.2818 Myr 



-200 -100 



100 200 300 400 

Radial distance (AU) 



100 200 



Fig. 4. Zoomed-in view on the surviving fragment in 
model 1 at four time instances. The color image shows the 
gas surface density in log g cm~^, the yellow curve outlines 
the fragment and the red circles marks the Hill radius. 



Section [2T3l The red circle represents the corresponding Hill 
radius of the fragment calculated as 



Ru = rf 



1 



Mf 



3 M, + Mf 



1/3 



(19) 



where is the stellar mass, rf is the orbital distance of 
the fragment and Mf is the mass of the fragment confined 
within the yellow curve. A mini-disk with a developed two- 
armed spiral structure can be seen around the fragment in 
the lower-right image. 

We used our fragment tracking algorithm to calculate 
physical properties of the surviving fragment. Figure [5] 
presents (a) the orbital distance of the fragment rf , (b) the 
mass of the fragment Mf and the mass confined within the 
Hill radius Mh (solid and dashed lines, respectively), (c) 
the radius of the fragment and the Hill radius Ru (soHd 
and dashed lines, respectively), and (d) the integrated grav- 
itational torque acting on the fragment T (solid line). The 
latter quantity is calculated as the sum of all individual 
torques r = — m(r, (j))d^/d(j) acting on the fragment, where 
m(r, (j)) is the gas mass in a cell with polar coordinates 
(r, 0) and ^ is the corresponding gravitational potential. 
The dotted line marks the zero torque. 

The orbital distance of the fragment varies in the 308- 
352 AU limits. The mean distance from the central star is 
ff = 330AU and eccentricity of the orbit is e ^ 0.07. Such 
low eccentricity orbits are typic al for companions formed 
by disk gra vitational instability (|Vorobvov fc Basul , l2010bl : 
lBossl . l2012D . The mass of the fragment varies in the 7-14 Mj 
limits. This wide scatter reflects either imperfections in the 
fragment tracking mechanism or continuous perturbations 
imposed onto the fragment by the circumfragment disk and 
spiral density wake (or both effects). These perturbations. 
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Time (Myr) 

2.272 2.274 2.276 2.278 2.280 



Time(Myr) 

2.272 2.274 2.276 2.278 2.21 



Hill radius 

fragment radius 



mass within Hill radius 
fragment mass 



2.270 2.272 2.274 2.276 2.278 

Time (Myr) 



2.274 2.276 2.278 2.280 

Time (Myr) 



Fig. 5. Characteristics of the surviving fragment in 
model 1: (a) orbital distance, (b) mass of the fragment 
(solid line) and that confined within the Hill radius (dashed 
line), (c) radius of the fragment (solid line) and the Hill ra- 
dius (dashed line), and (d) integrated gravitational torque 
acting on the fragment in units of 8.6 x 10^^ g cm^ s~^. 



however, do not lead to the fragment dispersal at least on 
time scales of our numerical simulations. The mean mass 
of the fragment calculated over two orbital periods shown 
in Figure [5] is Mf=ll Mj, which is still in the planetary 
mass regime. However, the mean mass confined within the 
Hill radius is Mh = 20.5 Mj, which implies that the frag- 
ment may accumulate some more material in the course of 
the evolution as it cools and contracts into planetary-sized 
object. On the other hand, if the fragment could not get 
rid of angular momentum, most of the material in the Hill 
sphere would ultimately land onto a circumfragment disk 
and the fragment may remain the planetary-mass regime. 

Figure [5t shows the radius_of the fragment and the Hill 
radius. The mean Hill radius Ru = 47 AU is greater than 
mean radius of the fragment Rf = 20 AU by more than 
a factor of 2. The mean scale height at the position of the 
planet is about Z = 40 AU, somewhat s mall er than the Hill 
radius . According to lCrida et al.l (|2006f ) and lKlev fc NelsonI 
(|2012f ), the gap opening criterion can be written as 



'iR^ 



50^ 



< 1.0, 



(20) 



where q = Mh/M* and Vt^ is the orbital frequency of the 
fragment. Substituting the corresponding mean values for 
i?H and Tf into equation (|2Q|) . noticing that q = 0.017 for 
= 1.2 M© and Mr = 20.5 Mj, and finally calculating 
u using equation (|8|) and mean disk temperature of 15 K 
at Tf = 330 AU, we estimated the left-hand side of equa- 
tion ([2Q|) to be ^ 0.9, thus marginally satisfying the gap 
opening criterion. The future orbital dynamics of the frag- 
ment can be predicted using the integrated gravitational 
torque acting on the fragment from the rest of the disk, T, 
shown in Figure [Hi. Evidently, the torque is mostly pos- 
itive, implying outward migration, and its mean value is 
T = 4.9 X 10""^ in units of 8.6 x 10"^° g cm^ s^^ 

We estimate the characteristic migration timescale us- 
ing the following simple analysis. A (small) change of the 
orbital distance dr^ of a fragment with mass Mf on a 



Keplerian orbit caused by a (small) change in the an- 
gular momentum of the fragment dL can be written as 
drf = 2 dL/MfVf^ where = {GM^/vfY^'^ is the angu- 
lar velocity of the fragment. The migration velocity of the 
fragment is then 



dVi 

~dt 



9 M± 

^ dt 



Noticing that dL/dt - 
can be calculated as 



(21) 



T, the characteristic migration time 



mg 



rf 



7mg 



L 



(22) 



Substituting the corresponding mean values for rf , Mf and 
T into equation ([22]) . and also noticing that = 1.2 Mq 
in model 1, we estimated the migration timescale to be of 
the order of 10 Myr. 

Since the fragment in model 1 opens a gap in the disk, a 
more appropriate estimate of the migration timescale may 
be that fflven by the v iscous diffusion time in the disk 
dLin fc PapaloizouL [19861 ) 
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(23) 



which yields essentially the same timescale of 10 Myr. The 
above estimates show that the fragment will remain on a 
wide orbit for a ti me period longer than the typical dis k 
lifetime of 2-3 Myr fetrom et al.[. [l989[ : [Haisch et al.[ . [2QQl[ ). 
implying that the fragment will finally turn into a massive 
GP or low-mass BD on a stable orbit of the order of 400 AU 
from the central object. 

5. Formation of an intermediate-mass brown dwarf 

In this section we describe the formation of a 43- Jupiter- 
mass BD around a 0.9 Mq star. Figure [6] shows a series of 
images of the gas surface density (logarithmic in g cm~^) 
in model 2 for the inner 1000 AU of our computational 
box. The time elapsed since the formation of the central 
protostar is indicated in each image. This model is charac- 
terized by initial gas surface density and angular velocity in 
the pre-stellar core that are independent of radial distance, 
i.e., S = const and Q = const. As a consequence, model 2 
has an elevated mass infall onto the disk, stronger grav- 
itational instability a nd more vi^o urous fragmentation in 
the Class phase (Vorobyov", '2012) than model 1, notwith- 
standing the fact that the prestellar core in the latter model 
is more massive. Gravitational instability in model 2, fueled 
by intense mass loading from the envelope, is so strong that 
the disk has broken into massive clumps linked with each 
other by long and dense filaments. During the course of 
the evolution, most fragments have migrated onto the star 
due to strong torques but one massive fragment manages 
to survive through the initial violent stage and settles onto 
a quasi-stable orbit at around 0.3 Myr. 

Figure [7] presents main characteristics of the surviv- 
ing fragment during approximately four orbital revolutions. 
The layout of the Figure is the same as that of Figure [H The 
orbital distance of the fragment varies in the 170-185 AU 
limits and the mean distance is ff = 178 AU. The fragment 
is characterized by a low eccentricity orbit, e ~ 0.04, some- 
what smaller than that of the fragment in model 1. The 
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Fig. 6. Gas surface density distribution in model 2 shown at 
various times since the formation of the central protostar. 
Only the inner 2000 x 2000 AU box is shown, the total 
computational region extends to 16000 AU. The scale bar 
is in log g cm~^. Note a fragment on a stable orbit in the 
bottom row. 
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Fig. 7. Same as Figure [5] only for model 2. 
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Fig. 8. Gas surface density distribution in model 3 shown at 
various times since the formation of the central protostar. 
Only the inner 2000 x 2000 AU box is shown, the total 
computational region extends to 20000 AU. The scale bar 
is in log g cm~^. Note two fragments on quasi-stable orbits 
in the third row of images. One of the fragments disperses at 
t ~ 1.55 Myr and the other survives to the end of numerical 
simulations {t = 1.8 Myr). 



The integrated gravitational torque acting on the frag- 
ment in model 2 is always positive and is somewhat stronger 
than that of model 1, possibly due to a higher mass of 
the former. The migration timescale calculated using equa- 
tion (|22]) is found to be of the order of 4 Myr, which is 
comparable to or even longer than the typical lifetime of 
the disk. We conclude that the fragment in model 2 has 
good chances to survive migration and settle on a wide or- 
bit evolving ultimately into an intermediate-mass BD. 



mean mass of the fragment is Mf = 43 Mj and the mean 
mass contained within the Hill radius is Mh = 60 Mj . This 
implies that the fragment will ultimately form an interme- 
diate mass BD, perhaps surrounded by its_own circum-BD 
disk. The mean radius of the fragment is = 21 AU and 
the mean Hill radius is Rn = 46 AU. The latter value is 
greater than the local scale height, Z ^ 24 AU, and the 
gap opening criterion (|2Q|) is satisfied. In general, the frag- 
ment in model 2 appears to be in a more perturbed state 
than that of model 1, perhaps, because the former is more 
massive, younger and less evolved. 



6. Attempted formation of multiple companions 

Notwithstanding the fact that protostellar disks in mod- 
els 1 and 2 exhibit multiple episodes of gravitational frag- 
mentation, only one fragment in each model has survived 
after 1.0 Myr of disk evolution. This raises the question 
of whether gravitational fragmentation can account for the 
formatio n of m ulti-companion systems similar to HR 8799 
(Marois et al.L [200 8. 2010), which has four planetary-mass 
objects on orbits at 15-70 AU from the central star. Below, 
we discuss this possibility. 
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Fig. 9. Main characteristics of fragment 1 (left column) 
and fragment 2 (right column) in model 3 before the least 
massive of them (fragment 2) has dispersed. In particu- 
lar, Panel a presents the orbital distance of the fragments, 
Panel b — masses of the fragments (solid lines) and those 
confined within the Hill radii (dashed lines), and Panel c — 
radii of the fragments (solid lines) and their Hill radii 
(dashed lines). 



Figure [8] presents a series of images of the gas surface 
density in model 3. The parameters of the model are listed 
in Table [1] and the time after the formation of the proto- 
star is indicated in each image. The initial evolution of the 
disk is characterized by vigourous fragmentation and sev- 
eral fragments at a time are usually present in the disk. 
By t = 1.1 Myr only two fragments survive and settle on 
quasi-stable orbits with only slightly different radial dis- 
tances from the star but with a 160° offset in azimuthal 
angle with respect to each other. However, after orbiting 
in unison for about 0.45 Myr, one of the fragments dis- 
perses at t ^ 1.55 Myr. The dispersed fragments leaves a 
a crescent-shaped density enhancement in the disk which 
can still be seen in Figure [8] at t = 1.6 Myr. The other 
fragment survives to the end of our numerical simulations 
(t = 1.8 Myr). 

Figure [9] shows main characteristics of the two surviv- 
ing fragments during a time period of 1.46-1.5 Myr, i.e., 
before one of the fragments has dispersed. In particular, 
the left-hand and right-hand panels belong to fragment 1 
(most massive) and fragment 2 (least massive), respectively. 
Panels (a) in Figure [9] present the orbital distance of the 
fragments, panels (b) — masses of the fragments (solid lines) 
and masses contained within the Hill radii of each fragment 
(dashed lines), and panels (c) — radii of the fragments (solid 
lines) and their Hill radii (dashed lines). 

Both fragments move on orbits that are less stable than 
those in models 1 and 2, perhaps due to continuing gravita- 
tional perturbation exerted on the fragments by spiral den- 
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Fig. 10. Same as Figure [5] only for model 3. 

sity wakes excited by both fragments in the disk. The mean 
orbital distances of fragment 1 and 2 are rf^i = 407 AU and 
ff^2 = 393 AU, respectively. The eccentricity of the orbits 
is also varying somewhat and the maximum eccentricity 
for fragments 1 and 2 are estimated to be approximately 
ei = 0.05 and 62 = 0.1, respectively. As a result, notable 
radial excursions are evident in the top panels of Figure [9l 
The masses of both fragments stay in the planetary-mass 
regime, though with significant variations reflecting their 
highly perturbed state, and the mean masses of fragments 1 
and 2 are Mf,i 5.9 Mj and Mf, 2=4.0 Mj, respectively. 
The mean masses contained within the JHill radii of frag- 
ment 1 and 2 are Mh,i=14.0 Mj and Mr, 2=8.9 Mj, re- 
spectively, indicating that if both fragments had survived 
they would have formed_massive GPs. The mean radii of 
fragments 1 and 2 are Rf^i=lS AU and Rf^2=l^ AU, re- 
spectively, and their mean Hill radii are i?H,i=48.5 AU and 
i?H,i=41 AU. Both values are greater than the local scale 
higiit Z ^ 25 AU. 

The likely reason why one of the fragments has dispersed 
is that this pair of fragments violate the criterion for orbital 
stab ility b etween two coplanar planets on circular orbits 
(Gladmad, [1993) 



Af > Af,, 



2V^i?H, 



M, 



(24) 



where Af = ff,2 — ^f,i = 14 AU and Rum is the mutual 
Hill radius defining the region in which gravitational force 
between two bodies is larger than the force on them due to 
the star 



^H,M = 



Mf,i + Mf,2 
3M* 



1/3 



TfA + rf ,2 



(25) 



Substituting the corresponding values into equations 
and (|25]) . one obtains Rum = 56 AU, Af^cr=193 AU, and 
Af <C Af^cr- Evidently, the orbits of the two fragments are 
unstable. Moreover, the mutual Hill radius is greater than 
the difference between the mean orbital distances of the 
fragments implying that the less massive fragment might 
not have withstood the disturbing tidal influence from the 
more massive counterpart. This effect might have been 
aided by insufficient numerical resolution of our logarithmic 
polar grid at a radial distance of the fragments (~ 400 AU). 
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We note that the value of Af^cr = 193 AU is of the order of 
the orbital distance for directly-imaged wide-orbit compan- 
ions (see forth column in Table [3]), which implies that the 
separation between companions in multicomponent systems 
should be comparable to their orbital distances. This makes 
the formation of such wide-orbit, multicomponent systems 
even more difficult. In any case, th e results of this numerical 
simulation and other studies (e.g. iBossl . I2Q1L) have demon- 
strated that gravitational fragmentation can account for the 
formation of multiple fragments at a time, but the question 
of whether these fragments can ultimately mature into a 
system with more than one companion in wide orbits is 
still open. 

After the less massive fragment dispersed at t ~ 
1.55 Myr, the other fragment settled on a quasi-stable orbit 
with mean radial distance f = 370 AU. The main charac- 
teristics of the surviving fragment are shown in Figure [TOl 
the layout of which is the same as that of Figure [H The 
last surviving fragment is characterized by the mean mass 
Mf = 11.0 Mj, which is similar to the total mass of the two 
fragments before one of them dispersed. This suggests that 
the surviving fragment has accreted most of the material re- 
leased by the destroyed fragment. Significant temporal vari- 
ations in the instantaneous mass of the surviving fragment 
are indicative of its highly pert urbed_st ate. The mean mass 
contained within the Hill radius is Mh = 19.0 Mq. The or- 
bit of the surviving fragment is characterized by rather low 
eccentricity, e = 0.02. The integrated torque acting on the 
fragment is positive and the estimated characteristic migra- 
tion time of the fragment is tmg = 6.9 Myr. We conclude 
that this fragment is likely to evolve into a massive GP or 
low-mass BD, depending on the amount of mass that will 
ultimately settle into a circumfragment disk. 

7. Characteristics of survived GP/BD embryos and 
comparison with observations 

We have run 60 models with the total integration time in 
each model exceeding 1.0 Myr after the formation of the 
central protostar. Protostellar disks in the majority of mod- 
els were sufficiently massive to experience vigorous gravita- 
tional fragmentation at radial distances greater than sev- 
eral tens of AU during the initial several hundred thousand 
years. The number of the fragments amounted to more than 
ten at a time. However, most of the fragments have ei- 
ther migrated through the inner computational boundary 
at 6 AU or got ejected from the computational domain into 
the intracluster medium or dispersed by tidal torques on 
time scales less than 1.0 Myr. Only six models out of 60 
revealed the survival of one of the fragments after 1.0 Myr 
of evolution. 

Figure [11] brings together six models that have demon- 
strated the formation of stable companions on wide orbits 
(to which we refer below as GP/BD embryos), showing 
for each model the gas surface density image (g cm~^ in 
log units) at the end of numerical simulations. The model 
number and time elapsed since the formation of the cen- 
tral protostar are indicated in each panel. Only the inner 
2000 X 2000 box is shown for each model. All six embryos 
possess their own circum-embryo disks, the masses of which 
are comparable to those of the parent embryos as implied 
by the mass contained within the Hill radius (see Table [2] 
below). In particular, circum-embryo disks in models 2 and 
5 exhibit a pronounces two- armed spiral structure. 



It is worth noting that the survived embryos and 
circum-embryo disks both exhibit a retrograde rotation 
with respect to that of the parental protostellar disk. This 
goes against expectations of havin g prograde rotation (e.g 
Bole v et aD . 120101 : IZhu et al.L [2012. ) borne out by the fact 
that the specific angular momentum of a fluid element in a 
rotationally supported disk increases with distance as r^-^. 
Our nume rical scheme perfo r ms w ell on standard numer- 
ical tests ([Vorobvov fc Basul . |2Q06[), including an an gular 
momentum conservation test (* Norman et all Il980[ ), and 
hence this phenomenon is unlikely to be caused by nu- 
merical reasons. However, signiflcant local deviations from 
Keplerian orbits ca n be present in st rongly gravitationally 
unstable disks (e.g. IVorobvovi |2010[ ). with a consequence 
that the speciflc angular momentum radial proflle may not 
locally reflect a rotationally supported disk. A visual anal- 
ysis of Figures [2l [6] and [8] demonstrates that both the 
retrograde and prograde fragments are present at the ini- 
tial stages of disk evolution, but the retrograde systems 
may be favoured for survival. For instance, the Type HI 
migration depends strongly on the flow pattern near the 
planet (|Peplihski et al.Ll2008[ ). The counterrotating circum- 
embryo disk may make the formation of horseshoe stream- 
lines near the corotation more diflicult, thus weakening the 
efficiency of Type HI migration. In any case, a more fo- 
cused study is needed to understand the phenomenon of 
retrograde rotation. 

Embryos in models 1-5 satisfy the gap opening crite- 
rion (|2Q|) and clear a well-defined gap in the protostellar 
disk. The embryo in model 6 is the least massive of all and 
the corresponding gap is less pronounced. We note that the 
gap is profoundly non-axisymmetric in models 2, 3, and 5, 
an effect that can in principle be used to infer the presence 
of massive GPs/BDs in the disk. 

We present the main characteristics of surviving em- 
bryos in Table [2j In particular, columns 1-11 list the 
model number, mass of the prestellar core Mc, ratio of ro- 
tational to gravitational energy in the prestellar core /3, 
mass of the protostar M*, mean mass of the embryo Mf, 
mean mass within the Hill radius Mh, mean orbital dis- 
tance of the embryo rf, orbital eccentricity e, mean ra- 
dius of the embryo i^f, mean Hill radius i?H, and migra- 
tion timescale tmg- For the purpose of comparison, we also 
provide in Table [3] the main characteristics of directly im- 
aged, wide-orbit (> 50 AU) companions to stars with mass 
0.08 M0 < M* < 2.1 Mq. We excluded lower-separation 
companions because we are interested only in GPs/BDs 
that can be formed by disk gravitational fragmentation, 
which is likely to occur at radial distances greater than 
50 AU. In particular, columns 2-5 list stellar masses M^^ 
masses of companions Mp, orbital distances of companions 
Tp, and stellar ages as compiled by the Extrasolar Planets 
Encyclopedia (http://exoplanet.eu). We have ordered the 
objects along the line of increasing stellar masses. 

The embryo masses in our models lie in the 3.5 — 43 Mj 
limits and cover the entire range of masses found for de- 
tected wide-orbit companions. There is no clear indica- 
tion that more massive cores tend to produce more mas- 
sive companions, though the range of companion- forming 
core masses considered in the present study (1.2-1.7 Mq) is 
quite narrow. There is a hint that models with higher P tend 
to produce more massive companions (e.g., models 2 and 5) 
and that BD embryos tend to orbit less massive protostars, 
but these tendencies need to be confirmed with a wider sam- 
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Fig. 11. Gas surface density distribution in models showing the formation of quasi-stable GP/BD embryos on wide 
orbits. The model number and time elapsed since the formation of the central protostar is indicated in each panel. Only 
the inner 2000 x 2000 AU box is shown. The scale bar is in log g cm~^. 

Table 2. Characteristics of survived embryos 



model Mc /3 M* Mf Mr ff 
(Mo) (%) (Mo) (Mj) (Mj) (AU) 



(AU) (AU) (Myr) 



1 


1.7 


0.56 


1.2 


11 


20.5 


330 


0.07 


20 


47 


10 


2 


1.2 


0.88 


0.9 


43 


60 


178 


0.04 


21 


46 


4 


3 


1.5 


0.56 


1.1 


11 


19 


370 


0.02 


20.5 


54 


6.9 


4 


1.4 


0.56 


1.0 


4.6 


11.3 


415 


0.03 


16.5 


46 


7 


5 


1.55 


1.27 


0.75 


27.5 


40.5 


180 


0.06 


19.5 


41 


3.9 


6 


1.4 


0.56 


0.95 


3.5 


5.0 


190 


0.05 


11.5 


20.5 


4.1 



pie of models. Our modelling suggests that the companion 
masses are likely to be determined by the disk/stellar prop- 
erties rather than by the properties of parental prestellar 
cores. The masses contained within the Hill radius are gen- 
erally a factor of 1.5-2.5 greater than those of the embryos. 
The final masses of the embryos are therefore expected to 
be somewhat higher, though the actual growth will be lim- 
ited by the angular mome ntum of the disk g as that the 
embryo is trying to accrete (jBolev et al.L 1201 Of ). 

The orbital distances of GP/BD embryos in our mod- 
els are confined in the 178-415 AU limits, whereas Table [3] 



indicates a wider range of orbital distances for directly im- 
aged companions, 15-1170 AU. While the upper limit on 
the orbital distance in our modeling (415 AU) may in- 
crease if we run higher-resolution numerical simulations 
(due to better resolution at large radial distances on the 
adopted logarithmically-spaced radial grid), the lower limit 
(178 AU) is not expected to change considerably given that 
all em bryos in our mode ls are in fact slowly migrating out- 
ward. 'Vorob vov fc Basul ([2010b) found GP embryos form- 
ing at distances of the order 50 AU but their models em- 
ployed a barotropic equation of state which is known to 
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Table 3. Characteristics of known wide-orbit GPs and BDs 



Star 




Mp 




Age 




(Mo) 


(Mj) 


(AU) 


(Myr) 



Oph 11 


0.16 


21 ± 3 


243 ± 55 


11 ± 2 


CHXR 73 


0.35 


-1 r'-l-S 

15tl 


210 


2 


DH Tau 


0.37 ±0.12 




330 


1 


CD-35 2722 


0.4 


31 ±8 


67 


100 


GSC 06214-00210 


0.6 


17 ±3 


320 


11 ± 2 


Ross 458(AB) 


0.6 


8.5 ±2.5 


1170 


475 ± 320 


GQ Lup 


0.7 


21.5 ±20.5 


103 


1 


IRXS J1609 


0.7 


^ 8 


330 


2-5 


CT Cha 


0.8 ±0.1 


17 


440 


2±1 


AB Pic 


0.8 


13.5 ±0.5 


260 


30 


HN Peg 


1.0 


16 ±9 


795 ± 15 


200 


HR 8799 


1.6 


5-10 


15-68 


20-150 


Fomalhaut 


2.1 


q+1.2 
"^-0.5 


119 


200 



facilitate disk fragmentation at small radial distances from 
the star. 

The dearth of embryos at orbital distances ^170 AU 
in our simulations cannot be explained by numerical res- 
olution effects and is intriguing as it conflicts with ob- 
servations. There are four objects in Table [3] that have 
companions orbiting the star at radial distances less than 
170 AU (CD-35 2722, GQ Lup, HR 8799, and Fomalhaut). 
If disk fragmentation cannot explain these objects, then 
they must have been formed by dynamical scattering of 
closely-packed companions leading to eje ction of the least 
mass i ve companio n onto a wide orbit (e.g. lScharf fc Menoul . 
120091 : IVeras et al.L 12009). The resulting orbits of these com- 
panions are expected to cover the full range of possible ec- 
centri cities, some of them may even become unbound with 
time (|Veras et al.L |2QQ9[ ). In addition, the probability to 
form wide-orbit companions via dynamical scattering de- 
clines quickl y with the increasing o rbital distance of the 
companions (|Scharf fc Menoul , l2QQ9f ). The eccentricities of 
our formed embryos lie in the 0.02-0.07 range and the only 
known eccentricity of a wide-orbit companion is that of 
Fomalhaut b, e = 0.11. Interestingly enough, the orbital 
distance of Fomalhaut b is 119 AU and fits into the scat- 
tering scenario for the formation of wide-orbit companions. 
The above analysis suggests a unified picture for the forma- 
tion of wide-orbit companions, in which objects at orbital 
distances from several tens to ~ 150 AU are preferentially 
formed by dynamical scattering and are characterized by 
a whole spectrum of eccentricities while those at greater 
orbital distances are mostly formed via disk fragmenta- 
tion and are characterized by low eccentricities. Definite 
measurement of eccentricities for other wide-orbit compan- 
ions should therefore clarify the formation mechanism of 
GPs/ BDs on wide orbits (see also iDodson- Robinson et aH 
120091 ). 

The minimum mass of a protostar that hosts a wide- 
orbit GP/BD embryo in our models is found to be = 



0.75 M0, whereas the directly-imaged wide-orbit compan- 
ions have host stars with masses extending down to the 
brown-dwarf-mass regime (e.g., UScoCTIO lOS. lBiar et aD . 

2008) . It appears that protostars with mass ^0.7 Mq pos- 
sess protostellar disks with mass that is insufficient to ex- 
perience gravitational fragmentation in the T Tauri phase. 
Fragmentation episodes in these disks are mostly confined 
to the embedded phase and are driven by mass infall from 
the parental core. As discussed in Section [3j fragments in 
such disks have little chance to escape fast inward migra- 
tion and are unhkely to form wide-orbit companions. We 
note that iBossI (|2Qllf ) reported the formation of compan- 
ions in disks around stars with mass < 0.7 M©, but the 
total integration time was limited to just 1000 yr and was 
insufficient to draw firm conclusions about the survivability 
of these companions. 

If disk fragmentation cannot explain wide-orbit com- 
panions around stars with mass ^0.7 M©, then a viable 
alternative is dynamical scattering of close-orbit compan- 
ions. However, Table [3] demonstrates that most of the stars 
with mass < 0.7 Mq have companions on orbital distances 
of the order of several hundred AU. Such large orbital dis- 
tances are difficult to explain in the framework of dynami- 
cal scattering since the number of scatter ed objects quickly 
declin es with increasing orbital distance (|Scharf fc Menom 

2009) . This inconsistency necessitates further research into 
this subject. 

Finally, we note that our modelling failed to produce 
systems with more than one stable companion. This is con- 
sistent with observations. The only known system with sev- 
eral wide-orbit companions is HR 8799 and even in this 
case only one companion has an orbital distance greater 
than 50 AU. The lack of multiple wide-orbit companions 
is likely caused by the fact that the orbital stability cri- 
terion (|24|) imposes strict limitations on the minimum or- 
bital distance between two stable companions Af^cr- For in- 
stance, in model 3 that showed an attempted formation of 
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two companions (but one of them finally dispersed), the cor- 
responding minimum distance was Af^cr = 193 AU, while 
the actual mean radial separation between the compan- 
ions was only 14 AU. Large values of Af^cr make the disk 
fragmentation scenario problematic for explaining multi- 
component systems as it may require very extended and 
hence massive disks. 



8. Conclusions 

We computed the gravitational collapse of prestellar cores 
with masses lying in the 0.1 Mq < < 1.8 Mq range and 
ratios of rotational to gravitational energy confined in the 
0.2% < p < 2.2% limits. The integration time in our nu- 
merical hydrodynamics simulations extended beyond 1 Myr 
after the formation of the central protostar and covered the 
entire embedded phase and part of the T Tauri phase of 
stellar evolution. We focused on models that showed disk 
gravitational fragmentation and, in particular, on models 
that revealed the formation of quasi-stable, giant planet 
(GP) and brown dwarf (BD) embryos on orbits greater 
than 50 AU (referred as wide-orbit companions). The typi- 
cal migration timescales of the embryos are comparable to 
or greater than the lifetime of a typical disk (2-3 Myr), 
which allows us to conclude that they will ultimately cool 
and contract into fully formed GPs or BDs. Insufficient nu- 
merical resolution does not allow us to resolve the formation 
of planetary-sized objects. 

While most of our models showed disk fragmentation, 
only 6 out of 60 models revealed the formation of wide- 
orbit companions. We compared the characteristics of our 
embryos with those of fully formed GPs and BDs obtained 
from direct imaging (http:/ /exoplanet.eu). The disk masses 
and radii provided below were calculated by time- averaging 
instantaneous values over the duration of the Class I phase, 
which is most relevant to studying disk fragmentation. Our 
findings can be summarized as follows. 

— Masses of wide-orbit companions lie in the 3.5 — 43 Mj 
limits and cover the entire range of masses found for 
directly-imaged GPs and BDs on orbits greater than 
50 AU. There is no clear indication that the mass of 
the companion depends on the mass and angular mo- 
mentum in the pre-stellar core, though a wider sample 
of companion-forming models is needed to draw firm 
conclusions. 

— The orbital distances of the companions found in our 
modelling lie in the 178-415 AU limits. This range of 
orbital distances is notably narrower than that found 
for directly-imaged companions, a few AU-1170 AU. 
While the upper found limit (415 AU) may increase if we 
consider higher resolution simulation^ the lower limit 
(178 AU) cannot be explained by resolution effects and 
is intriguing as it conflicts the observations. We propose 
that companions at orbital separations from a few tens 
to 150 AU are likely to form via dynamical scattering of 
closely-packed companions onto wide orbits, while com- 
panions at larger orbital distances are predominantly 
formed via disk fragmentation. Definite measurement 
of eccentricities as a function of orbital distance should 
clarify the formation mechanism of GPs/BDs on wide 



^ Numerical resolution deteriorates at large distances on our 
logarithmically spaced radial grid. 



orbits, because the disk fragmentation scenario tends to 
produce companions on low eccentricities, e < 0.1. 

— Our numerical simulations did not produce multiple 
companions on wide orbits. Although model 3 revealed 
an attempted formation of two companions at orbital 
distances ~ 190 — 210 AU, one of the fragments dis- 
persed in less than 0.5 Myr. The likely reason for this 
failure is that the radial separation between the two 
companions, ^14 AU, was much smaller than the re- 
q uired 190 AU acco rding to the orbital stability criterion 
of iGladmanI (|l993f ). The formation of multi-component 
systems would require very extended and hence massive 
disks, which are statistically rare. The only known sys- 
tem with several wide-orbit companions, HR 8799, has 
only one companion at a distance greater than 50 AU 
and is likely to form via dynamical scattering. 

— The minimum mass of a companion-hosting star found 
in our modelling is 0.75 M©, whereas the correspond- 
ing value for directly imaged systems extends down to 
the BD-mass regime. It is likely that disk fragmentation 
in systems with stellar mass < 0.7 Mq is mostly con- 
fined to the embedded phase of star formation, in which 
mechanisms causing migration/ejection/destruction of 
the fragments are strong and the likelihood for survival 
of the fragments is low. 

— Disk gravitational fragmentation does not automati- 
cally guarantee the formation of wide-orbit compan- 
ions. Most of the fragments do not stay on wide or- 
bits for more than several orbital periods. The major- 
ity are torqued into the disk inner region and through 
the sink cell (6 AU) by gravitational torques from spiral 
arms or other fragments, a few fragments may be ejected 
from the disk via many-bo d y grav it ational interaction 
(Sta matellos fc WhitworThl . l2009al : iBasu fc Vorobvovl 
2012), and some fragments may be t idally dispersed at 
various radial distances in the disk (iBolev et al.L 201Qt 
iNavakshinL 120101 : IVorobvovL l2011al : IZhu et al.L I2012D . 
The fragments that pass through the sink cell may be 
comple tely destroyed and accreted onto the formi ng star 
(IVorobvov Basul . I2006L l2010al: |Machida eliD, 1201 ll) 
or lose their gaseous envelopes and form terrestrial cores 
or icy giants if th e dust sed i mentation timesca l e was suf- 
ficiently short (fNavakshinl . l2010l : iBolev et al.l . |2010[ ). 

— The minimum disk mass at which fragmentation can 

fr 

take place is found to be = 0.07 Mq for models 

with (3 ^ 1.2%. The value of may increase by al- 
most a factor of 2 for disks formed from pre-stellar cores 
with lower angular momentum, 0.2% < P < 1.2%. Our 

fr 

found values of are in reasonable agreement with 
previous est imates of the critical disk mass, ^0.1 Mq 
(jRice et al.L I2003L iMaver et all I2007D . 

— The minimum disk mass that is required for the forma- 

c.f . 

tion of wide-orbit companions is found to be = 
0.21 M0, a factor of 3 greater than the minimum disk 

mass for gravitational fragmentation = 0.07 Mq. 
The minimum disk radius that is required to produce 
wide-orbit companions is 370 AU, more than a factor of 
two greater than that required for disk fragmentation 
to occur. 

The overall low probability of the fragment sur- 
vival and the efficient migration/destruction mechanisms 
found in the present study both support the recently 
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proposed disruption/ d ownsizing hypothesis (|Bolev et al.L 
I2Q1QI : [Nayakshinl, |2Q1Q[ ) . Our res ults are also in agreement 
with the recent study of IZhu et~ aL (2012), who studied nu- 
merically gravitational fragmentation in two-dimensional 
disks and found that only 3 out of 13 fragments escaped fast 
inward migration. In the present paper, we have ignored a 
possible motion of the central star due to a combined grav- 
itational potential of the non-axisymmtric disk. Analytic 
studies suggest that the stellar wobbling may increase the 
strength of the odd s piral modes, and in particular that of 
the m = 1 mode (e.g. lShu et al.L[l99Qf ). which may increase 
the gravitational torque exerted on the fragments and make 
it even more difficult for the fragments to survive. On the 
ot her hand, focused numeri cal hydrodynamic simulations 
bv lMichael fc DurisenI (|2Q1Q| ) found no evidence for the en- 
hanced m = 1 mode in a 0.14 Mq disk around a solar- mass 
star. The effect of stellar wobbling in more massive disks, 
relevant for the present results, needs to be considered in a 
future study. 
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